General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



N83>11011 


(ti&SA>TH-eJ983) OU TB£ 0BI6XU Of 
BJiLATIV13TIC PABTICLBS AMD GAMMA-BAYS IB 
UOASAES (MASA) 18 p HC A02/MP A01 CSCL 03B 

Unclas 
G3/93 35286 


IVV\SA 

Technical Memorandum 83983 


On the Origin of Relativistic 
Particles and Gamma-Rays in 
Quasars 



J. Protheroe and D. Kazanas 


AUGUST 1982 


National Aeronautics and 
Space Administration 

Goddard Space Flight Center 

Greenbelt. Maryland 20771 



ON THE ORIGIN OF RELATIVISTIC PARTICLES AND GAMMA-RAYS IN QUASARS 


R. J. Protheroe^ and D. Kazanas^ 

Laboratory For High Energy Astrophysics 
NASA/Goddard Space Flight Center 

Re- :elved 1982 May 19 

Abstract 

A model Tor a class of quasars and active galactic nuclei Is proposed In which 
a shock around a massive black hole randomizes the Infall kinetic energy of spheri- 
cally accreting matter producing a non-thermal spectrum of high energy protons. It 
is suggested that these protons are responsible for the secondary production (via w* 
decay) of the radio emitting high energy electrons and also of high energy Y-rays 
(via decay and inverse Compton interactions of the electrons). The correlation 
between radio and Y-ray emission Implied by the model is in good agreement with 
observations of 3C273. Observation of the flux of high energy neutrinos from quasars 
may provide a test for the model. 


To be published in the Astrophvsica 1 Journal. 


Ulso University of Maryland 


2 


I. Introduction 

Accretion onto compact objects has long been considered as an energy source for 

many astrophysical phenomena, both on stellar and galactic scales. Due to Its 

potentially high efficiency (possibly up to 10%), It was one of the first mechanisms 

considered to power quasars (Lyndon Bell 1969) and Is thought to be the only 

mechanism capable of satisfying their enormous power requironents. However, It was 

realized early on that matter falling freely Into a black hole need not radiate all 

Its potential energy (Shapiro 1973). Accretion on to a black hole can be efficient 

only If the directed Infall motion can be randomized near the black hole, e.g., by 

dissipation In an accretion disk (e.g. Shields, 1978), or by dissipation of magnetic 

fields and relativistic particles convected with the Infalling plasma (e.g. 

Meszaros, 1975). Radio, X-ray and Y-ray observations show that some quasars emit 

roughly the same energy per decade (i.e. per logarithmic energy Interval, d log^QE) 

from radio to X-ray energies and. In the case of 3C273, to y-ray energies. Since 

Infall In a gravitational potential is unlikely to release more than ^ 0.1 mc^ per 

particle, not sufficient to produce the observed y-rays, we suggest that the 

relativistic particles required to account for the observed radio, X and y-ray 

emission may result from acceleration in shocks via a first-order acceleration 

(up to 98t; Axford 1981) 

mechanism. This mechanism, efflclenty^ln transforming directed kinetic energy Into 
relativistic particle energy, has already been proposed to account for the observed 
cosmic ray spectrum (Bell 1978a, b; Axford, Lear and Skadron, 1977; Blandford and 
Ostrlker, 1978; Cowsik and Lee, 1981). For strong shocks, this process naturally 
produces a power law spectrum with index -2, i.e., equal energy per decade. 

The nuclear interaction of protons accelerated by this process results In 
production of y-rays and secondary electrons responsible for synchrotron and Inverse 
Cjmpton emission. X-ray emission will be due to the Interaction of these relatl- 
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vistic particles with the ambient medium, photons, and magnetic f1 Ids. In this 
paper we will concentrate solely on the production of relativist.,, particles and 
r-rays in quasars. 

II. The Model 

The model employs a shock around a massive black hole which serves both to 
randomize the infall kinetic energy of the spherically accreting gas and also to pro- 
duce a non-thermal distribution of high energy (relativistic) protons. A thermal 
distribution o^ particles may also be produced at the shock and hence radiate by thermal 
breiTKstrahlung; however we consider here only the radiation resulting from the non- 
thermal distribution which may actually contain most of the energy (Axford 1981). 

Because of the high photon energy densities (5^ 200 erg cm" for 3C273 if the source 
region has dimensions less 

than 10^^ cm) and the magnetic fields present, electrons will suffer severe Compton 
and synchrotron losses, effectively prohibiting their direct acceleration. Direct 
production of electron-positron pairs by interactions of protons with photons may 
become important at very high energies (Blumenthal, 1970) but we shall consider here 
only those secondary electtons resulting from nuclear collisions of the non-thermal 

proton distribution (via pion production) as opposed to those produced by 
thermal distributions of protons (Marscher, Vestrand and Scott, 1980). Gamma- 
ray production via inverse Compton scattering of the relativistic electrons 
(Jones, l'^79) and via 71® decay are considered. 

The question of the existence of a shock around a black hole is important 
and deserves further discussion. It has been argued, on general grounds, that 
accretion shocks can exist only if the gravitating object has a solid surface 
to provide the necessary boundary condition requiring the deceleration of the 
supersonically accreted matter. A black hole, since it does not provide such 
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a surface, apparently violates this condition. The situation, however, is 
different if relativistic particles (protons, in the present case) of sufficient 
energy density, c, are present around the black hole. These particles having 
kinetic energies much higher than their gravitational potential do not readily 
fall into the black hole. Moreover outward magnetic field gradients will tend 
to nirror then, thus further inhibiting their accretion. Thus the relativistic 
particles can provide sufficienc pressure to satisfy the Rankine-Hugoniot con- 
ditions across the shock. The confinement of the particles is then achieved by 
the ram pressure of the accreting material, i.e. «. “i; pv^. The position of the shock 
is then fixed by this condition and by the balance betw»:n injection and removal (by 
nuclear collisions) of the relativistic particles (protons). The initial formation 
of the shock can be attributed either to convection of ancient relativistic particles 
or to dissipation of magnetic fields as suggested in the literature (McCrea 1956; 
Scharlemann 1931). 

The spectrum of protons accelerated by shock has a power law momentum 
dependence, i.e., Qp(p) = Kpp"^ (protons GeV-^ c s-M- For a strong shock 
(compression ratio = 4), r * 2. The relativistic protons are assumed to be trapped 
within 3 radius comparable to the shock radius by magnetic fields which are tied to 
the infalling matter. If the relativistic protons were at confined to the central 
region of the quasar, spallation of nuclei in the outer gas clouds responsible for 
the line emission in quasars would result in an overabundance of boron (Baldwin et 
al., 1977), contrary to observation. In addition, the matter density in this region 
is assumed to be sufficiently high that the relativistic protons are depleted mainly 
by nuclear interactions. This condition sets a lower limit on the accretion rate 
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for a given black hole mass which however is not very restrictive. Nuclear 
interactions throughout the volume within the shock radius will result in a spectrum 
of pions at production given by: 

Q,(E) “ ^ |E’,E)dE' (pions GeV-i s“M (1) 

where i(i^(E',E) is the probability of a proton of energy E* producing a pion of 
energy E in a single nuclear interaction and is obtained from accelerator data. The 
factor r/(r-l) takes into account multiple interactions of the primary protons 
(Protheroe, 1981). 

The production spectrum of y-rays, resulting from ir® decay (Stecker 

1971) has been calculated and is plotted in Figure 1(a) where it is normalized to 
the proton production spectrum at high energies, KpE"^. The production spectrum of 
secondary electrons Qg(E) (e'‘‘e") from it-y-e decay has been calculated and is also 
plotted in Figure 1(a). Above 1 GeV, Qg(E) » Q^^(E) - 0.’5 Qp(E) for r = 2; the 

remaining energy goes into neutrino production, Q^(E) » 0.7 Qp(E). To obtain the 
absolute normalization of Qp(E) for 3C273, we shall consider its observed radio- 
optical spectrum. Based on the survey of flux measurements on 3C273 by Ulrich 
(1981) and using Hq = 60 km s"^ Mpc“^, we find that the luminosity emitted per decade i 
approximately constant at ' 1.5 x 10^® erg s"^ from 3 x 10^^ to 10 ^^ and the total 
luminosity is about 2 x 10^^ erg s"^ integrated over all wavelengths. If we attribute 
the flat (F^ ' v'M part of the spectrum from the radio to the visible to synchrotron 
radiation, then we can relate the sychrotron luminosity to the total electron 
production rate: 


L 


synch 


rt ? -1 ^ 

«^nl0 E (E) (erg s ^ decade"'"). 


( 2 ) 
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where a Is the ratio of energy densities In the magnetic and radiation fields. 

Qg^°^(E) Is the rate of production of electrons at high energies (electrons erg"^ s"M 
and Is equal to Q^(E) plus the production spectrum of electrons and positrons resulting 
from photon-photon pair production Interactions. 

Assuming the electrons are tripped within the shock radius, R, and lose energy by 
inverse Compton and synchrotron radiation, the ambient spectrum of electrons, Ng(E), Is 
determined by 

Oe“V) • 0. (3) 

.here Q,“*(E) ■ ♦ 7 |Q,^(E') * Q,c(E')) ♦^(E') dE'. (4) 

Here, dE/dt - - 3.97 x lO-^ (1 + a) e»*g whe»*^* Is the radiation energy 

density, taken to be Ljq^j^/hR^c. The second term of equation (4) Is the energy 
spectrum at production of electrons from pair production and Involves the y-ray pro- 
duction spectrum resulting from inverse Compton Interactions as well as that from ir° 

decay. ♦ (E) is the probability that a y-ray of energy E will produce an electron- 

yy 

positron pair by annihilating with a lower energy photon before escaping from the 
source (discussed In detail later). For Inverse Compton Interactions, the y-ray 
production rate Is, 


Qlc(E) =£/“ 2 o(E,c,E’)cde}dE' (5) 

E/4y ^ 

where y * E'/mc^, Is the number density of photons In the radiation field 

(photons cm"^ erg~M, and o(E,e,E') Is the differential cross section for an Inverse 
Compton int jractlon in which an electron of energy E' Interacts with a photon of energy 


e, boosting its energy to E. 
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We have solved equation (3) for the two extreme cases where the source Is 


optically thin • 0) or thick *1) to y-rey pair production Interactions. 
At high energies* but below the Kleln-Nishina limit, we obtain, for r ■ 2; 


0,“'(E) ■« 


0,IE) • 0) 

. Q.^(E)| . 1) 


( 6 ) 


The resulting Inverse Compton y-ray production spectra are plotted In Figure 1(») for 
these two cases and a = 1. (Not using the Kleln-Nishina cross section may Introduce an 
error of about a factor of 2 at 10 GeV.) 

The observed y-ray spectrum will be the total y-ray production spectrum 
(Q + Q^p) modified by photon-photon pair production Interactions. To calculate this 
effect, we must know the radius of the object, R. We shall obtain R from the frequency 
below which the synchrotron radiation becomes optically thick to synchrotron self 
absorption. The number density of relativistic electrons of energy E within the source 
Is 

00 

n^(E) = {(IPTHF (electrons erg"^ cm"^) . (7) 


Using Qg^®^(E) from equation (2), and substituting for dE/dt, we obtain 
nglE) - n^ E*^, where n^ = (4.9 x lO^ Lj^^g^)/(oLyg^j^R). For this spectrum, the 
frequency below which synchrotron self absorption becomes Important, v^. Is given by 
Pachoiczyk (1970), 


2Cj(SCg) 


2/7 


n 2/7 q 5/7 


( 8 ) 


where Is the component of magnetic field perpendicular to tne line of sight 

(B^ » .9<B>), S • j R is the average distance through the source, Cj* 6.3 X 10^®, and 

Cg = 7.8 X 10"^^. The synchrotron self absorption does not result In a sharp spectral 
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break In 3C273. This may be understood In terms of non-uniform radial distributions of 
electron density and magnetic field (Condon and Dressel, 1973). For 3C273, appears 
to be In the range 3x10^^ - 10^^ Hz, hence we find - 5x10^^ - 10^^ cm and 

- 14 - 75 gauss. 

We can now use these estimates of the radius to obtain the photon density In the 
X-ray region from the observed X-ray flux data and hence the effects of pair production 
Interactions on the r-ray spectrum. The optical depth to v-rays of energy E Is 

t(E) “ S / ''pjj(c) 0 (E»e) de (9) 

where Op^le) Is the number density of photons of energy e, and o(E,e) Is the diffe- 
rentia! cross Section for pair production Interactions of Y-<*ays of energy E with 
photons of energy e, (Heltler, 1954). We have used X-ray data of Worrall et al. (1979) 

and Primlnl et al. (1979) and the assumed radius to obtain n^. (e). 

ph 

The energy spectnjm of y-rays escaping from the source Is then down by a 
factor, (1 - e”^)/T, from that at production. For 3C273 this factor Is plotted as a 
function of energy In Figure 2(a) for various assumed radii. We see that the source Is 
optically thick to y-rays above ' 100 MeV for radii In the range obtained earlier. We 
have therefore plotted In figure l^b) the y-ray spectrum at production (Q + Q,«) 

TTY 

for * 1 and for values of o ranging from 1 to 1/16. The observed y-ray luminosity 
(Bignami et al., 1981; Hermsen et al., 1981) Is compared In Figure 2(b) with our 
prediction taking Into account photon-photon absorption and Is fouro to be In good 
agreement with the prediction for a - 0.1 and R - lO^^ cm. This radius Is consistent 
with that we obtained earlier by considering the radio spectrum some 12 decades lower 
in energy. 
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III. Discussion 

We have presented above the main features and predictions of a model for a class 

of quasars and active galactic nuclei, gauged to 3C273. The model, deductively 

Introduced, ^ocuses on the high energy particles In these objects \vh1ch are responsible 

for the radio and y-ny emission. It predicts a correlation between these two extreme 

parts of the spectrum which appears to be In good agreement with observations of 3C273, 

the only quasar In which both radio emission and high energy yrdys have been 

observed. The model provides estimates of the central source radius and Its magnetic 

47 -1 

field. If the total observed luminosity = 2x10 erg s ) Is die to 

gravitational energy release at the surface of the shock, and the relativistic particle 
pressure balances the accretion ram pressure pv^, then the radius Inferred above 
Implies a mass of 10^° Mg, and an accretion rate of 'lO^ Mg yr~^. 

The estimated radius Is compatible with variability time scales of the order of 
several months which have been observed {Worrall et al. 1979). Variability reported on 
much shorter time scales {- 1 day) (Marshall, Warwick and Pounds, 1981) would clearly 
be Incompatible with our model. However, the statistical significance of these data Is 
low and they may also be compatible with a much longer Lime scale variability* In any 
case, variability on time scales much less than R/c has been observed in the 
quasar NRAO 140 (Marscher and Broderick, 1981) and has been attributed to relativistic 
beamlmj effects. It Is Important to point out here that a radius nf - 1 light day 
would render 3C273 completely opaque to v-rays energies greater than a few MeV, 
contrary to observation. However, emission of X-rays and v-rays from the jet (McBreen, 
1979; Bassani and Dean, 1981, Blandford and Konigl, 1979; see also Konigl, 1981) has 
been suggested and might be a possible remedy If the observed X-rays are beamed with 
the Y-rays. 

In the present model, pair produc^ion by high energy y-rays and production of 
electrons by pi on decay throughout the source volum'- avoids the Compton Catastrophe 
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(Hoyle, Burbldge and Sargent, 1966) and allot^s a higher energy density In the radiation' 
field than In the iragnetlc field (I.e. a < 1). The model also provides a situation 
similar to that considered In the literature for the formation of jets (Smith et a1. 
1980), I.e., a hot gas (the relativistic r<>*t1c1es) underneath a cold gas (the 
Infalling cold gas). 

A prediction of the model would be the emission of high energy v-rays frois other 
quasars and active galactic nuclei. Unfortunate', sources of luminosity similar to 
3C273 are typically -lO timei farther away and difficult to detect In v-rays at 
present. The present model may also apply to Seyfert galaxies and. If so, many of 
these objects may be detectable In v-rays (Bassani end Dean, 1981). Howevrr, their 
radii are at present uncertain and do not preclude their being optically tnick to y- 
rays. 

Finally, the model predicts a higher flux of high energy neutrinos than In models 
where electrons are directly accelerated. We predict a flux of 1.4 x 10"^® x (l+2o) 
neutrinos (all types) cm“^ s"^ above 1 TeV for 3C273. This may be within the detection 
capability of the proposed DUMAND neutrino telescope (Stenger, 1981) and thus may 
provide a test of the modtl. 

In conc’.;si*>n, we have presented a model ^or a class of quasars and active 
galactic nuclei which naturally provides the relativistic particles needed to account 
for their spect**a. It also avoids the Compton catastrophe, Inherent In most models, 
and predicts a correlation between the spectrum at radio and y-ray energies. We have 
successfully applied this model to 3C273. 

We are Indebted to R. F. Mushotzky for Introducing us to extragalactic astronoiqy 
and for providing advice and criticism which has been valuable In understanding some 
of the problems associated with extragalactic objects. We would also like to thank 
0. C. Ellison, R. Ramaty, R. A. Shafer, R. Sllberberg, W. A. Stein, C. M. Urry, and A. 
W. Wolfendale for stimulating discussions. 
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Figure Captions 

Figure 1. (a) Production spectra of electrons (Qg)i If-rays from decay 

rays fum Inverse Compton Interactions for a » 1 (Qxc^ relative to the 
production rate of protons* Separate curves are given for corresponding 
to • 0 and 1. (b) Total production rate of ^or • 1 relative 

to the radio synchrotron luminosity for different values of a* 

Figure 2 * (a) Reduction factor applied to predicted optically thin y-ray spectra for 

3C273 for various assumed radii in order to obtain the expected y-ray 
emission, (b) Gamma-ray luminosity of 3C273 based on the observed spectrum 
(Hermsen et al., 1981'> compared with the luminosity predicted for a * 1 and 
a * I/IO ind for the range of radii obtained as described in the text. 
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